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SUMMARY 


, . i e ^y namic ’ P ro Pulsion, and mass models for a generic, horizontal-takeoff, single-stage- 
to-orhn (SSTO) configuration are presented which are suitable for use in point mass as well as 
batch and real-time six degree-of-freedom simulations. The simulations can be used to investigate 
ascent performance issues and to allow research, refinement, and evaluation of integrated guidance/ 
ilightypropulsion/thermal control systems, design concepts, and methodologies for SSTO 
missions. Aerodynamic force and moment coefficients are given as functions of angle of attack 
Mach number, and control surface deflections. The model data were estimated by using a 
subsonic/supersonic panel code and a hypersonic local surface inclination code. Thrust coefficient 
and engine specific impulse were estimated using a two-dimensional forebody, inlet, nozzle code 
and a one-dimensional combustor code and are given as functions of Mach number, dynamic 
pressure, and fuel equivalence ratio. Rigid-body mass moments of inertia and center of gravity 
location are functions of vehicle weight which is in turn a function of fuel flow. 

INTRODUCTION 

The basic objective of the National Aero-Space Plane (NASP) Program is the development 
of a manned honzonta 1 takeoff and landing, single-stage-to-orbit (SSTO), airbreathing launch 
vehicle Realization of this concept is theoretically possible through advanced technology in the 
areas of aerodynamics, materials, structures, flight control, and propulsion (ref. 1). Critical to the 
development of the vehicle is the integration of these disciplines to arrive at a technology mix 
which maximizes propulsive efficiency while minimizing aerodynamic heating and structural loads. 



The studies reported in references 1-3 for a conical accelerator configuration, tested 
experimentally from low subsonic to high hypersonic Mach numbers (refs. 4-6), show that system 
integration involving aerodynamics, propulsion, and guidance and control early in me design can 
significantly reduce the fuel fraction to orbit. The mathematical models developed for these 
investigations contained sufficient detail to allow study of trim drag reduction concepts, guidance 
and control strategies, and vehicle performance issues. In reference 2 the models were modified 
slightly to allow thrust vectoring and active center-of-gravity (c.g.) location control to be 
investigated as ways to reduce excessive fuel consumption caused by drag due to elevons used to 
trim the vehicle. 

The purpose of this report is to present the aerodynamic, propulsion, and mass 
mathematical models used in the investigations reported in references 1*3. Aerodynamic 
longitudinal and lateral-directional force and moment coefficients are given as functions of angle of 
attack; Mach number, and elevon, rudder, and canard deflections. These model data were 
estimated using a subsonic/supersonic panel code and a hypersonic local surface inclination code. 
Thrust coefficients and specific impulse are given as functions of Mach number, dynamic pressure, 
and fuel equivalence ratio, and were estimated using a two-dimensional forebody, inlet, and nozzle 
code and a one-dimensional combustor code. The moments of inertia and center of gravity 
position are given as functions of vehicle weight which in turn varies with fuel flow. 

SYMBOLS AND ABBREVIATIONS 

b lateral-directional reference length, span, ft 

c longitudinal reference length, mean aerodynamic chord, ft 


Gn, total drag coefficient, , n.d. 

D CJ S re | 


drag increment 

Cp. drag increment coefficient for basic vehicle, 


C D,da 

C D,de 

C D,dr 

C D,dc 


drag increment A 

drag increment coefficient for right elevon, =rr , n.d 

q ^rcf 


drag increment 


drag increment coefficient for left elevon, , n.d 

q ^rcf 


drag increment coefficient for rudder, Q , n.d. 

Q Or 


drag increment 

q S rc f 


drag increment coefficient for canard, — » n.d. 


dra g increment 
Q Sref 



total lift coefficient, 


__L_ 

q S r ei 


n.d. 


2 



C La 

^L,da 

C L,de 

C, . 
L,dc 


lift increment coefficient for basic vehicle, increment n ^ 

q S re f 


lift increment coefficient for right elevon, iift increment n ^ 

q Srcf 


lift increment coefficient for left elevon, li ft increment ^ n j 

q Srcf 


lift increment coefficient 


for canard. 



n.d. 


total side force coefficient, — ~ , n.d. 

Q ^rcl 



side force with sideslip derivative for basic vehicle 


dCy 

"air 


n.d. 


C 

i,da 


side force increment coefficient for right elevon, force increment n ^ 

q Srcf 


c 

u Y,de 


side force increment coefficient for left elevon Side force increment n j 

q S r cf ’ ' ’ 



side force increment coefficient for rudder, Si de force increment n j 

qS re f ’ ‘ * 


C 


1 


total rolling moment coefficient , — ^ , n.d. 

q S r cfb 


c 


1 


p 


rolling moment with 


sideslip derivative for basic vehicle. 


dCi 

dP * 


n.d. 



rolling moment increment coefficient for right elevon ffl.pment increment d 

q Srcfb 



rolling moment increment coefficient for left elevon D10me D t increment , 

q S re fb ’ n a 
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l,dr 


rolling moment increment coefficient for rudder, mom ^y n crement ^ n ^ 


aq 


rolling moment with roll rate dynamic derivative, p ^ , n.d. 


2 V 


3 C 

rolling moment with yaw rate dynamic derivative, ^ n.d. 

5 2V 


M 


C total pitching moment coefficient, =r-r , n.d. 

m F q S re fC 


ma 


pitching moment increment coefficient for basic vehicle, moni SD (increment ^ n j 

q o rc fC 


'm,da 


pitching moment increment coefficient for right elevon, increment ^ n d 

q ^refC 


'm,de 


pitching moment increment coefficient for left elevon, increment , n ^ 

q ^refC 


'm,dr 


pitching moment increment coefficient for rudder, m t >m ent incremen t ^ n j 

q ^rcfC 


"m,dc 


pitching moment increment coefficient for canard, increment ^ n d. 

q ^rcfC 


m 


a C m 

pitching moment with pitch rate dynamic derivative, — t— , n.d 

a-^ 

2 V 


N 


total yawing moment coefficient , =-= — — , n.d. 

q S rc fb 


’P 


dCn 

yawing moment with sideslip derivative for basic vehicle, , n.d. 

ap 
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c , 

n,da 

C , 
n,de 

C n,dr 


yawing moment increment coefficient for right elevon, i noment increment n d 

q S ref b 


yawing moment increment coefficient for left elevon increment j 

q Srefb ’ na 


yawing moment increment coefficient for rudder moment increment n d 

q S re fb 




c T 

Isp 

L 

M 

N 

Minrc 
L, D, Y 
X, Y, Z 

a 

P 

<t> 

M 

q 


yawing moment with roll rate dynamic derivative. 


aCn 
2 V 


n.d. 


yawing moment with yaw rate dynamic derivative, n d 

r b ’ 

2 


thrust coefficient, - -^ st , ft^ 

engine specific impulse, sec 
aerodynamic rolling moment at center of gravity, ft-lb 
aerodynamic pitching moment at center of gravity, ft-lb 
aerodynamic yawing moment at center of gravity, ft-lb 
aerodynamic pitching moment at moment reference center, ft-lb 
total aerodynamic lift, drag, and side forces respectively, lb 
total aerodynamic x, y, and z body axes forces respectively, lb 
angle of attack, deg 

sideslip angle, rad 

engine fuel equivalence ratio, n.d. 

Mach number, n.d. 

dynamic pressure, lb/ft 
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2 

S f reference area, theoretical wing area, ft 

T engine net thrust, lb 

V vehicle freestream velocity, ft/sec 

W fuel flow rate, lb/sec 

W vehicle weight, lb 

Wo initial value of vehicle weight, lb 

Wcon weight of fuel consumed, lb 

longitudinal distance from moment reference center to vehicle c.g., positive aft, ft 

2 

I X x, Iyy> Izz roll, pitch, and yaw moments of inertia respectively, slg-ft 

APAS Aerodynamic Preliminary Analysis System 

c.g. Vehicle center of gravity 

G&C Guidance and Control 

HABP Hypersonic Arbitrary Body Program 

MACH Mach number, n.d. 

mrc wind tunnel moment reference center 

n.d. nondimensional 

POST Program to Optimize Simulated Trajectories 

SSTO Single-Stage-To-Orbit 

UDP Unified Distributed Panel program 

VEHICLE DESCRIPTION 

The planform and side view of the study vehicle are shown in figure 1 , and the geometric 
characteristics are given in table I. This aircraft is a full scale version of the wind tunnel test 
configuration reported in reference 4. A sizing analysis discussed below yielded a full-scale gross 
weight of 300,000 lbs and an overall fuselage length of 200 ft. 

The fuselage-centerline-mounted wing has conventional, independently controllable, 
trailing edge elevons with their hinge line perpendicular to the fuselage centerline. Deflections of 
the elevons are measured with respect to the hinge line, and positive deflections are with the trailing 
edge up. The fuselage is modelled as an axisymmetric conical forebody, a cylindrical engine 
nacelle section, and a cone frustrum engine nozzle section. A fuselage-centerline-mounted vertical 
tail has a full span rudder with its hinge line at 25 percent chord from the trailing edge. Deflections 
of the rudder are measured with respect to its hinge line, and positive deflections are with the 
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ttailing edge left. The small canards have a 6 percent thick symmetrical 65A series airfoil and are 
deployed only at subsonic speeds. Deflections of the canard are measured relative to the fuselage 
centerline, and positive deflections are with the trailing edge down. 


AERODYNAMICS MODEL 


A subsomc/supersomc/hypersonic analysis code developed jointly by NASA Langley and 
Rockwell International Inc., referred to as the Aerodynamic Preliminary Analysis System (APAS) 
was used to predict the longitudinal-and lateral-directional force and moment coefficients of the 
study configuration. The APAS program is described conceptually in Appendix A. 

estimated the basic vehicle lift and drag increment coefficients and the side force 
sideslip derivative as functions of angle of attack and Mach number. The program estimated the 
lilt, drag, and side force increment coefficients for the right and left elevons, the rudder, and the 
canard as functions of angle of attack, surface deflection, and Mach number. Changes in lift drag 
and side force due to body angular rates and aerodynamic coupling between control surface ' " ' 

deflections and sideslip are assumed negligible. 


APAS estimated the basic vehicle roll and yaw moment sideslip derivatives, the pitch 
moment increment coefficient, and the roll, pitch and yaw dynamic derivatives (for roll, pitch and 
yaw rates) as functions of augle of attack and Mach number. The roll and yaw moment increment 
coefficients for right and left elevon and rudder and the pitch increment coefficients for right and 

et * xf V u 0, rud u der ’ ^ d canard were estimated as functions of angle of attack, surface deflections 
and Mach number. These quantities are given relative to the moment reference center. The total 
moments relative to the c.g. are obtained by adding those caused by lift, drag, and side force to the 

ab ?\f t j uantltie ; s - Aerodynamic moment coupling between control surface deflections sideslip 
and body angular rates is assumed negligible. 


Data were computed at Mach numbers of 0.3, 0.7, 0.9, 1.5, 2.5, 4.0, 6.0, 10.0, 15.0, 
20.0, and 24.2 and angles of attack of -1.0°, 0.0°, 1.0°, 2.0°, 4.0°, 6.0°, 8.0°, 10.0°, and 12.0 C 
At each Mach number and angle of attack combination, coefficients were generated for a range of 
deflections of the right elevon, the left elevon, and the rudder, each taken separately. Deflections 

of -20.0 , -10.0°, 0.0°, 10.0° and 20.0° were used for each surface. At the three subsonic Mach 
numbers and nine angles of attack above, data were estimated for canard deflections of -10.0°, 

-5.0° 0.0° 5.0° and 10.0°. The APAS Unified Distributed Panel (UDP) code was used at Mach 
numbers of 0.3, 0.7, 0.9 and 1.5, and the Hypersonic Arbitrary Body Program (HABP) code 
was used at Mach 2.5 and above. The results of the computations are shown in figures 2-32. 


Model Implementation 

i • ^ m P^ ementat ^ n °f aerodynamic model in a computer simulation requires a numerical 
algorithm to extract the increment coefficients and derivatives. The drag, lift, and side force and 

rolling, pitching, and yawing moment increment coefficients and derivatives for the basic vehicle 
are determined as functions of angle of attack and Mach number. The increment coefficients 
caused by control surface deflections are determined as functions of angle of attack, surface 
deflections, and Mach number, and are added to the basic vehicle increments to form the total 
aerodynamic force and moment coefficient. These coefficients are used to compute the forces and 
moments relative to the moment reference center. These moments are then corrected to the c g 
1 he following sections discuss the computation of these quantities 
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Dra g Force Computation 

The basic vehicle drag increment coefficient, is shown in figure 2 as a function of 
angle of attack and Mach number. The drag increment coefficients for the right elevon, d , left 
elevon, C D d , rudder, C D df , and canard, C D dc , are shown in figures 3-6, respectively, as 

functions of angle of attack, deflection angle, and Mach number. These coefficients are summed to 
obtain the total drag coefficient as 

C D = C Da + C D,de + C D,da + C D,dr + C D,dc 

and the drag is given by 


D_qS ref C D 


Lift Force Computation 


The basic vehicle lift increment coefficient, C La’ 


is shown in figure 7 as a function of 


angle of attack and Mach number. The increment coefficients for the right elevon, d , left 

elevon, dg , and canard, d{ ., are shown in figures 8-10, respectively, as functions of angle 

of attack, deflection angle, and Mach number. Changes in lift due to rudder deflections are 
neglected. These coefficients are summed to obtain the total coefficient as 


C L “ C La + C L,de + C L,da + C L,dc 


The lift is then given by 


L “ q S ref C L 


Side Force Computation 

The basic vehicle side force due to sideslip derivative, C v , is shown in figure 1 1 as a 

Y P 

function of angle of attack and Mach number. The increment coefficients for the right elevon, 

Cy da , left elevon, Cy d , and rudder, Cy ^ , are shown in figures 12-14, respectively, as 

functions of angle of attack, deflection angle, and Mach number. The total side force coefficient is 
given by 


C Y ~ C Yp P + C Y,da + C Y,de + C Y,dr 
where the side slip angle (3 is in radians. The side force is then given by 


Y_qS ref C Y 
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Rolling Moment Computation 


The basic vehicle rolling moment with sideslip derivative, C, , is shown in figure 15 as a 

P 

function of angle of attack and Mach number. The rolling moment increment coefficients for right 
elevon, Cj^ a , left elevon, Cj ^ and rudder, Cj are shown in figures 16-18 respectively as 

functions of angle of attack, deflection angle, and Mach number. The rolling moment dynamic 
derivatives for roll rate, C, , and yaw rate, C, , are shown in figures 19 and 20, respectively, as 

P r 

functions of angle of attack and Mach number. The total rolling moment coefficient is given by 


C 1 ~ C lp P + C l,da + C l,de 


+ C l,dr + C 1 



where P is in radians, and the terms (y^-) and (^y) are the computed nondimensional roll and 
yaw rates. The rolling moment is then given by 


L-qbS ref C l 

For this configuration there is no rolling moment caused by the z-axis force or the y-axis force 
because the c.g. and the moment reference center lie on the x-axis. 


Pitching Moment Computation 

The pitching moment increment coefficient relative to the moment reference center for the 
basic vehicle, C ma> is shown in figure 21 as a function of angle of attack and Mach number. The 

pitching moment increment coefficients for the right elevon, C- . , left elevon, C J rudder 
m,dr’ anc ^ canan *> ^m,dc ’ m s ^ own * n figures 22-25, respectively, as functions of angle of 
attack, deflection angle, and Mach number. The pitch rate dynamic derivative, C m , is shown in 

figure 26 as a function of angle of attack and Mach number. The pitching moment coefficient 
relative to the moment reference center is given by 



- r +c +c 
ma in, da rn,de 


+ C 


m,dr 


+ C . +C 
m,dc m 



The term (2 y ) ^e computed nondimensional pitch rate. The pitching moment relative to the 
moment reference center is given by 

M =q c S r C 
mrc 4 ref m 

For the study configuration the aerodynamic pitching moment about the c.g. is equal to the sum of 
the pitching moment relative to the moment reference center plus the pitching moment caused by 
the z-axis force acting through the x-distance from the c.g. to the moment reference center as 
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where the z-axis force is given by 




Z = - D sin a - L cos a 

For the study configuration there is no pitching moment caused by the x-axis force because the c.g. 
and the moment reference center lie on the x-axis. 

Yawing Moment Computation 

The yawing moment with sideslip derivative relative to the moment reference center for the 
basic vehicle, C , is shown in figure 27 as a function of angle of attack and Mach number. The 

"p 

increment coefficients for the right elevon, C n ^ left elevon, C n ^ and rudder, C n c j f , are 

shown in figures 28-30, respectively, as functions of angle of attack, deflection angle, and Mach 
number. The yawing moment dynamic derivatives for roll rate, C n , and yaw rate, C n , are 

P r 

shown in figures 31 and 32, respectively, as functions of angle of attack and Mach number. The 
yawing moment coefficient relative to the moment reference center is given by 


C n = C n p P + C „,da + C n,de + C „,dr + C „ p (I?) + C n r (fv) 
and the aerodynamic yawing moment relative to the moment reference center is given by 


N =qbS r C 
mrc H ref n 

For the study configuration the total aerodynamic yawing moment about the c.g. is equal to the 
sum of the yawing moment relative to the moment reference center plus the yawing moment caused 
by the side force acting through the x-distance from the c.g. to the moment reference center as 


N = N + x Y 
mrc eg 


where the y-axis aerodynamic force is given above. For this configuration there is no yawing 
moment caused by the x-axis force because the c.g. and the moment reference center lie on the 
x-axis. 


MASS MODEL 

The vehicle mass model is based on the assumption of rigid structure and time varying 
mass, center of gravity, and moments of inertia. The total mass of the vehicle, its c.g. location, 
and its moments of inertia vary as fuel is consumed. Fuel slosh is not considered, and the 
products of inertia are assumed negligible. 
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Model Generation 


The vehicle takeoff weight and size were determined based on an iterative method (refs. 
1,7). This method involves assuming vehicle takeoff mass, vehicle size, propulsion system, and 
aerodynamics and running the three degree-of-freedom Progr am to Optimize Simulated Trajectories 
(POST) simulation (ref. 8) to orbit to obtain a final mass fraction. The mass fraction is then used 
as input to an empirically based vehicle sizing program which outputs new values for mass and 
size. This process is repeated until no significant change is output from the sizing program. 

The sizing analysis yielded a vehicle gross weight of 300,000 lbs and an overall fuselage 
length of 200 feet. The c.g. location and the moments of inertia were estimated at ten values of 
vehicle weight ranging from 295,000 lb to 140,000 lb which corresponded to takeoff and engine 
shut down conditions, respectively. The moments of inertia are consistent with values given in 
reference 9 for several low aspect ratio military and research aircraft and the Concorde supersonic 
transport. The precomputed c.g. location, x , measured positive aft from the wind tunnel 

moment reference center, and the vehicle moments of inertia I YY , I YY , and I™ are shown in 
figure 33 as functions of vehicle weight. 

Model Implementation 

As in the aerodynamic model, implementation of the mass model in a computer simulation 
requires an algorithm to extract the mass data. The weight of fuel consumed, Wcon , is obtained 
by integrating the fuel flow rate over the time, t t, the engines are thrusting as 



The vehicle weight is then given by 


W - Wq - Wcon 


where Wo is the initial weight of the vehicle. 

The c.g. position, x cg , and the moments of inertia are then determined as a function of the 
vehicle weight. It is assumed that the c.g. moves only along the body x-axis as fuel is consumed. 

PROPULSION MODEL 

The propulsion model for the study configuration was developed using a two-dimensional 
forebody, inlet, and nozzle code with a one-dimensional combustor code. The code analyzes a 
vehicle nose-to-tail with the forebody, inlet, and nozzle flows computed by a two-dimensional 
perfect gas code. The combustor performance is computed assuming one-dimensional flow with 
hydrogen combustion, utilizing the equilibrium chemisuy combustion model of reference 10 The 
combustor code is based on an update of the cycle analysis method described in reference 1 1. The 
forebody, inlet, and nozzle forces are determined using a two-dimensional finite difference code 
based on the shock fitting method of reference 12. The boundary layer effects are determined by 
the methods described in references 13 and 14. 

The thrust coefficient (thrust per unit dynamic pressure) and specific impulse were 
estimated at specified values of Mach number, dynamic pressure, and equivalence ratio. The thrust 
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as determined from this model is net thrust, and it acts along the body x-axis. The effects of angle 
of attack and sideslip, body angular rates, and control surface deflections on thrust and specific 
impulse are considered negligible for the the present configuration. 

Thrust coefficients were computed at Mach numbers of 0.0, 0.3, 0.5, 0.7, 0.9, 0.95, 1.0, 

1.5, 2.0, 3.0, 3.5, 4.0, 6.0, 8.0, 10.0, 15.0, 20.0, and 25.0; five dynamic pressures ranging 

from 0 to 5000 lb/ft 2 ; and five equivalence ratios ranging from 0.0 to 10.0. Similarly, specific 
impulse values were estimated at Mach numbers of 0.0, 0.3, 0.5, 0.7, 0.9, 1.0, 1.5, 2.0, 3.0, 

3.5, 4.0, 6.0, 8.0, 10.0, 15.0, 20.0, and 25.0 ; five dynamic pressures ranging from 0 to 

5000 lb/ft 2 ; and six equivalence ratios ranging from 0.0 to 100.0. The results of the computations 
are shown in figures 34 and 35 respectively. 

Model Implementation 

As above, implementation of the propulsion model in a computer simulation requires an 
algorithm to extract the precomputed data. Using a linear interpolation/extrapolation scheme, the 
thrust coefficient and specific impulse are determined as functions of fuel equivalence ratio, 
dynamic pressure, and Mach number. The following sections discuss the computation of the 
thrust, specific impulse, and fuel flow rate. 

Thrust Computation and Throttling 

The thrust coefficient, Gp as shown in figure 34, is a function of equivalence ratio, 
dynamic pressure, and Mach number. The thrust is calculated as 

T = q Gp 

Engine throttling is simulated by varying the equivalence ratio. Equivalence ratio of unity 
corresponds to maximum fuel efficiency, and values greater than unity give more thrust but use 
disproportionately more fuel. 

Specific Impulse and Fuel Flow Computation 

The specific impulse, I sp , is shown in figure 35 as a function of Mach number, dynamic 
pressure, and equivalence ratio. Once the thrust is computed and the specific impulse is 
determined the fuel flow rate is computed as 



and is integrated to determine the weight of fuel consumed, Wcon • 

MODEL VALIDATION AND VERIFICATION 

The aerodynamic model presented is based on established conceptual design estimation 
methods. Good comparisons with available experimental aerodynamic data have been obtained 
(for example, in ref. 4, 5, 6, and 15.) 

The APAS program predictions compare well with Space Shuttle orbiter data book values 
as shown by C. I. Cruz at the Langley Research Center. There do not exist experimental data on 
control effectiveness for the study configuration; however, the APAS estimates of control 
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effectiveness have been validated in reference 15 for the Space Shuttle orbiter, the X-15, and a 
hypersonic research airplane (ref. 16). These comparisons indicate that APAS estimates of control 
effectiveness are generally satisfactory for conceptual design studies. 

The codes used to generate the propulsion model have been validated independently by 
comparisons with estimates from computational fluid dynamics codes and with limited 
experimental data. The sharp conical forebody and axisymmetric inlet, combustor, and nozzle of 
the study configuration improve the accuracy of the assumptions used in developing the model. 

For example, the cross flow around the forebody at angles of attack and side slip tends to keep the 
difference in thrust small between the windward and leeward sides of the vehicle. 

The vehicle mass characteristics and size were determined based on an established iterative 
method that involved the trajectory, the aerodynamics, and the propulsion system (refs. 1, 7). The 
sizing program uses empirical data, and the vehicle moments of inertia were made a function of 
vehicle weight based on the trends given in reference 9. 

The models presented have been used extensively in single-stage-to-orbit, three degree -of- 
freedom and six degree-of-freedom simulations (refs. 2, 3). In these cases active guidance and 
control systems were simulated and vehicle performance issues were studied. In reference 2 the 
model was modified slightly to allow thrust vectoring and active c.g. location control to be 
investigated as ways to reduce excessive fuel consumption caused by drag due to elevons used to 
trim the vehicle. In all cases the models gave reasonable and consistent results. 

Based on the above considerations, the models presented are suitable for use in point mass 
as well as batch and pilot-in-the-loop, six degree-of-freedom, horizontal takeoff and landing, 
ascent-to-orbit analyses and simulations. The models can be used to investigate trajectory 
optimization, guidance algorithms, stability augmentation system issues, drag reduction concepts, 
and to allow research, refinement, and evaluation of integrated guidance/flight/propulsion/thermal 
control system concepts and design methodologies for horizontal-takeoff single-stage to orbit 
missions. 


CONCLUDING REMARKS 

Aerodynamic, propulsion, and mass mathematical models for a generic, horizontal takeoff 
and landing, single-stage-to-orbit, conical accelerator configuration are presented. The models can 
be used to represent the vehicle in point mass as well as batch and pilot-in-the-loop, six 
degree-of-freedom, horizontal takeoff and landing, ascent-to-orbit analyses and simulations. 

These models can be used to investigate trajectory optimization, guidance algorithms, stability 
augmentation system issues, drag reduction concepts, handling qualities, and to allow research, 
refinement, and evaluation of integrated guidance/flight/propulsion/thermal control system concepts 
and design methodologies for horizontal-takeoff, single-stage to orbit missions. Aerodynamic and 
propulsion forces and moments are given as functions of local flow conditions, aerodynamic 
control surface deflections, body rates, and propulsion parameters. Mass properties are given in 
terms of vehicle weight which in turn varies with fuel flow. The models are based on accepted 
conceptual design methods and have been validated and verified separately and together in 
simulations. 
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APPENDIX A 

Aerodynamic Preliminary Analysis System (APAS) 

The FORTRAN-based Aerodynamic Preliminary Analysis System (APAS) (ref. 17)utilizes 
the Unified Distributed Panel (UDP) program for analysis in the subsonic/supersonic speed range 
and an enhanced version of the Hypersonic Arbitrary Body Program (HABP) in the high 
supersonic and hypersonic speed range. APAS is contained in an interactive shell which allows 
the user to quickly generate, modify, and analyze arbitrary vehicle geometries. 

The UDP program (ref. 18) estimates fuselage contributions to vehicle aerodynamics using 
slender body theory. The wing and tail contributions are evaluated using a system of vortex panels 
(non linear vortex effects are estimated using the Polhamus suction analogy.) The components of 
viscous drag, wave drag, and base drag are predicted by a combination of theoretical and empirical 
methods in UDP. The implementation of the vortex panel theory in UDP requires streamwise 
panel edges; thus, engineering approximations had to be made when defining the geometry at the 
wing/body and tail/body junctures, particularly in the aft region of the vehicle where the wing and 
tail join the truncated cone section of the fuselage. 

The HABP module (ref. 19) of APAS computes inviscid and viscous solutions on each 
finite element of the APAS model and integrates the pressure and shear forces to obtain total force 
and moment coefficients for the configuration. Pressures on the impact surfaces (surfaces whose 
normal vector has a component pointing into the flow) of the APAS model were approximated 
using the tangent-cone (fuselage component) and tangent-wedge (wing and tail components) 
methods. Prandtl-Myer expansion was assumed for all shadow surfaces. The viscous shear 
forces were estimated using the reference enthalpy method. 
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TABLE I.- GEOMETRIC CHARACTERISTICS OF CONFIGURATION 


Wing: 

Reference area (includes area projected to fuselage centerline), ft 2 3603 

Aspect ratio I oo 

S P an - ft 60.0 

Leading edge sweep angle, deg 75.97 

Trailing edge sweep angle, deg 0 0 

Mean aerodynamic chord, ft g() () 

Airfoil section diamond 

Airfoil thickness to chord ratio, percent 4 0 

Incidence angle, deg 0 0 

Dihedral angle, deg 0 o 

Wing flap (clevon): 

Area each, ft 2 92 3 

Chord (constant), ft 7 22 

Inboard section span location, ft 15 0 

Outboard section span location, ft 27 78 

Vertical tail, body centerline: 

Exposed area, ft 2 ^45 7 

Theoretical area, ft 2 1248 8 

32.48 

Leading edge sweep angle, deg 70 0 

Trailing edge sweep angle, deg 38 13 

Airfoil section diamond 

Airfoil thickness to chord ratio, percent 4 0 
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TABLE I. CONCLUDED 


Rudder: 


Area, ft^ 

Span, ft 

Chord to vertical tail chord ratio, percent 

Canard: 

Exposed area, 

Theoretical aspect ratio 

Span, ft 

Leading edge sweep angle, deg 

Trailing edge sweep angle, deg 

Airfoil section 

Incidence angle, deg 

Dihedral angle, deg 


161.4 

22.8 

25.0 


154.3 

5.48 

33.6 

16.0 

0.0 

NACA 65A006 

0.0 

0.0 


Axisymmetric fuselage: 

Theoretical length, ft 

Cone half angle, deg 

Cylinder radius (maximum), ft 

Cylinder length, ft 

Boattail half angle, deg 

Boattail length, ft 

Moment reference center, ft 


200.0 
. 5.0 

12.87 

12.88 
9.0 

40.0 

24.01 
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Figure 3.- Drag increment coefficient for right elevon 
attack, deflection angle, and Mach number. 
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Figure 3.- Concluded. 
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Figure 17.- Rolling moment increment coefficient for left elevon as a 
function of angle of attack, deflection angle, and Mach number. 
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Figure 18.- Rolling moment increment coefficient for rudder as a 
function of angle of attack, deflection angle, and Mach number 
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Figure 18.- Concluded. 
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Figure 19.- Rolling moment with roll rate dynamic derivative for basic 
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Figure 22.- Pitching moment increment coefficient for right elevon as a 
function of angle of attack, deflection angle, and Mach number. 
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Figure 23.- Pitching moment increment coefficient for left elevon as a 
function of angle of attack, deflection angle, and Mach number. 
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Figure 24.- Continued. 
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Figure 25.- Pitching moment increment coefficient for canard as 
function of angle of attack, deflection angle, and Mach number. 
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Figure 26.- Pitching moment with pitch rate dynamic derivative as a 
function of angle of attack and Mach number. 
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Figure 27.- Yawing moment with sideslip derivative for basic vehicle 
as a function of angle of attack and Mach number. 
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Figure 28.- Yawing moment increment coefficient for right elevon as a 
function of angle of attack, deflection angle, and Mach number. 
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Figure 28.- Continued. 
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Figure 28.- Continued. 
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Figure 28.- Concluded. 



Deflection 
angle, deg 

— □ — -20 


0 2 4 6 8 10 12 lc 

Angle of attack, a, deg 
(a) Mach number = 0.3 

Figure 29.- Yawing moment increment coefficient for left elevon as a 
function of angle of attack, deflection angle, and Mach number. 
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Figure 29.- Continued. 
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Figure 32.- Yawing moment with yaw rate dynamic derivative as a 
function of angle of attack and Mach number. 



Figure 33.- Moments of inertia and center of gravity (relative 
moment reference center) as a function of vehicle weight. 
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Figure 34.- Continued. 
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(o) Mach number = 10.0 



(p) Mach number = 15.0 


Figure 34.- Continued. 
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(q) Mach number = 20.0 



(r) Mach number = 25.0 
Figure 34.- Concluded. 
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(a) Mach number = 0.0 



Fuel equivalence ratio, <J> 
(b) Mach number = 0.3 


Figure 35.- Engine specific impulse as a function of fuel equivalence 
ratio, dynamic pressure, and Mach number. 
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(e) Mach number = 0.9 



(f) Mach number = 1 .0 


Figure 35.- Continued. 
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Fuel equivalence ratio, <|> 
(h) Mach number = 2.0 


Figure 35.- Continued. 
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Figure 35.- Concluded. 
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